Polar Behavior in a Magnetic Perovskite Via A-Site Size Disorder 
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We elucidate a mechanism for obtaining polar behavior in magnetic perovskites based on A-site 
disorder and demonstrate this mechanism by density functional calculations for the double perovskite 
(La,Lu)MnNi06 with Lu concentrations at and below 50%. We show that this material combines 
polar behavior and ferromagnetism. The mechanism is quite general and may be applicable to a 
wide range of magnetic perovskites. 

PACS numbers: 77.84.Dy,75.50.Dd 



There is great interest in materials that combine mag- 
netism and polar behavior, especially multiferroics with 
both ferromagnetism and ferroelectricityfii^i^ii^i^i^ While 
there are many fcrromagncts and many ferroclectrics, 
there remarkably few materials combining the two prop- 
erties. Arguments have been made explaining this ap- 
parent incompatibility in the case of perovskite ABO3, 
oxidesi^ In essence, it is because the best perovskite 
magnets have magnetic ions on the B site, while fer- 
roclectrics, such as BaTiOa, usually have B site ions 
with no d electrons. Here we elucidate a mccahanism 
for inducing polar behavior in magnetic perovskites and 
demonstrate this mechanism by density functional cal- 
culations for a ferromagnetic Ni-Mn double perovskite. 
This mechanism is quite general and may be applicable 
to a wide range of magnetic perovskites. 

Perovskite lattice instabilities are often understood us- 
ing the tolerance factor t = (ro -t- rA)/V^{r -\- tb), 
where ro, ta, and tb are the O, A-site and i?-site ionic 
radii, respectively^^ Ferroclectrics, such as BaTiOs and 
KNbOa, have t > 1, indicating the that B site ion is too 
small for its site in the ideal cubic structure. In the ferro- 
electric ground state of these so-called _B-sitc driven ma- 
terials, this ion off-centers, aided by hybridization with O 
states. There is another important class of ferroelectric 
perovskites, so called A-site driven materials. In these, t 
is normally less than unity, and the ferroelectricity is from 
off-ccntcring of A-site ions. This family includes the tech- 
nologically important Pb based piczoclcctrics and relaxor 
ferroclectrics. The essential physics is lone pair stere- 
ochemistry, specifically hybridization of Pb and Bi 6p 
states with O p statesi^ This class includes the few known 
magnetic ferroclectrics with strong ferroelectric proper- 
ties, e.g. BiFeOa, BiMnOg, and PbYOgi^iiaiiiii^iH With- 
out Pb or Bi, t < 1 perovskite structures generally derive 
from BOq octahedral tilts and not A-sitc off-ccntcring. 

Ions with d electrons are generally larger than d'^ ions. 
The majority of magnetic perovskites have t < I, with 
lattice structures based on tilts of the BOq octahedra and 
not ferroelectricity. However, first principles calculations 
have shown that, while these materials have tilted ground 



states, if the octahedra are prevented from tilting, strong 
ferroelectricity may result, with an energy intermedi- 
ate between the ideal cubic perovskite structure and the 
ground state structure, but closer to the later fiiii^ii^iii 
The role of Pb and Bi is then to shift the balance be- 
tween these states to yield ferroelectricity. Here we use 
a different approach to shift the balance between these 
states based on A-site size disorder J^ii^iii This is ap- 
plied rare-earth double perovskites, i?2MnNi06 where 
we obtain polar behavior combined with ferromagnetism 
for mixtures of large and small rare earth ions. The 
motivation for this choice is that the charge difference 
(5(5=2 between Mn^+ and Ni^+ and their size difference 
(Shannon radii^^ r]y[n-i+=0.67A, rNi2+=0.83A) indicates 
S-site ordering into the double perovskite structurcyi^ 
and that La2MnNi06 and Bi2MnNi06 are known to 
form and to be ferromagneticj ^°'^^i^^'^'^i^^'^^ These two 
compounds were previously studied by first principles 
calculations 

We used the local density approximation (LDA) in 
the general potential linearized augmented planewave 
(LAPW) method)^ with well converged basis sets in- 
cluding local orbitalsi^ The LAPW sphere radii were 2.0 
ao for La and Lu, 1.9 oq for Ni and Mn and 1.55 oq for 
Oj^ Mn and Ni atoms were placed in supercells with the 
double perovskite (rock-salt) ordering, and various or- 
derings of the A-site ions. The primary results reported 
here were done with 40 atom supercells. However, the 
lattice parameter was determined by relaxation of a 10 
atom cell of compositions LaLuMnNiOg- No symmetry 
was imposed, either for the 10 atom or 40 atom cells, but 
for the 40 atom cells the lattice parameters were held 
fixed at their pseudocubic values as determined from the 
relaxation of the 10 atom ferromagnetic cell, which as 
shown in Fig. [T] was 3.75 A. For both the parallel and 
antiparallel spin alignments, relaxation of the 10 atom 
cells yielded structures with off-centering of the Lu ions, 
i.e. polar structures, even though perovskite R25 type 
tilts are allowed in this cell. 

Our electronic structure for (La,Lu)MnNi06 is sim- 
ilar to those previously found for the Bi and La 
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FIG. 1: Energy of 10 atom LaLuMnNiOg cells vs. pseudocu- 
bic lattice parameter (a/2) for parallel (FM) and antiparallel 
(AF) Ni and Mn moments. Note that the FM alignment is 
strongly favored. 
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FIG. 3: (color online) Schematic depiction of the electronic 
structure of LaLuMnNiOe showing the dominant superex- 
change coupling by the arrow. 
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FIG. 2: (color online) Electronic density of states of a 10 atom 
ferromagnetic LaLuMnNiOe cell (a=3.77A). The peaks at ~ 
-3.7 eV and ~ 3.7 eV are the Lu and La / states, respectively. 
The projections are onto LAPW spheres. 



analogues and show high spin Mn'''+ and Ni^+. The 
LDA density of states for the relaxed ferromagnetic struc- 
ture at a=3.77A is shown in Fig. [2l and schematically 
in Fig. [31 For the various supercells we find either very 
small gaps or small band overlaps in the LDA depending 
on the exact crystal structure. We also performed some 
LDA-j-U calculations (not shown). These yield insulat- 
ing band structures, with gaps depending on U. Near 
metallicity is highly unfavorable for ferroelectricity, as 
it means strong electronic dielectric screening that will 
weaken the Coulomb interactions. Nonetheless, we use 
the LDA to avoid the ambiguity associated with an ad- 
justable parameter (U) and expect that the prediction of 
polar behavior will be robust. 

The ferromagnetism is due to the fact that with par- 
allel alignments of the Mn and Ni spins there is a strong 



TABLE I: Lu displacements for the various supercells. | < 
5 > I is the magnitude of the average Lu displacement, and 
< I (5 1 > is the average magnitude. 
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cross-gap hybridization of the unoccupied eg majority 
states of Mn"*"*" with the occupied majority eg states 
of Ni^+, leading to a ferromagnetic coupling consis- 
tent also with the Goodenough-Kanamori rules i'^^'^^i'^^ 
This ferromagnetic superexchange is particularly strong 
because in perovskites the strongest coupling through 
the near linear B-O-B bonds is via eg - Pa hopping. 
The strength of this coupling is evident from the siz- 
able crystal field splittings''^ of both the Mn and Ni d 
states (Fig. [2]). Thus the hybridization between oc- 
cupied and unoccupied eg orbitals, allowed for ferro- 
magnetic alignment, but not for antiferromagnetic align- 
ment, strongly favors ferromagnetism. Since the hopping 
is mediated by O this is not direct exchange, like the 
weak ferromagnetic coupling generally associated with 
the Goodenough-Kanamori ferromagnetism, but is a con- 
ventional strong superexchangei^2i2^i^ 

While our relaxed structure for the 10 atom cell is po- 
lar, this size supercell favors ferroelectricity (ferroelec- 
tricity arises from a zone center instability, while the 
competing tilt modes occur at the zone boundary and 
a 10 atom cell restricts tilts to the i?-point). We there- 
fore performed supercell calculations with a 2 x 2 x 2 40 
atom supercell. This cell is doubled along the [001], [Oil] 
and [111] directions and so allows arbitrary mixtures of 
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FIG. 4: (color online) Cation displacements with respect to 
their O cages in a 40 atom relaxed supercells of composition 
(Lao.5Luo.5)MnNi06. 



tilt instabilities and accommodates the observed Glazer 
patterns of perovskite tilt systemsi^ 

As mentioned, all supercells were for the double per- 
ovskite structure, i.e. rock-salt ordering of the i3-site Ni 
and Mn ions. For the A-site, we considered four arrange- 
ments of the Lu and La at a 50-50 composition as well 
as one supercell each for 5/8-3/8 and 3/4-1/4 composi- 
tions. The specific cells at the 50-50 composition were 
(1) rock-salt ordering of La and Lu ("G" in the follow- 
ing), (2) (001) layers of La and Lu ("A"), (3) lines of 
La and Lu along [001] and ordering c(2 x 2) in plane 
("C") and (4) a cell maximizing like near neighbors (La 
at corner and edge centers, Lu at face and body centers, 
denoted "P"). The 5/8-3/8 supercell was constructed by 
replacing one Lu by La in the cell "G" , while the 75-25 
supercell was made by substitution of two Lu by La in 
the same "G" supercell. 

Fig. 3] shows the cation positions in the lowest energy 
relaxed structures for the 50% Lu supercells with respect 
to the centers of their nearest neighbor O cage (nearest 
12 O atoms for the A-site atoms, and nearest 6 O for the 
_B-sites ions). The average displacements of the Lu for 
the various cells are sumarized in Table IH As may be 
seen, in all cases Lu strongly off-centers, by on average 
^ O.4A for the 50-50 supercells and slightly less for the 
lower concentration cells. Interestingly, unlike other A- 
site driven perovskite ferroelectrics;^ there is very little 
off-centering of the B-site ions. In fact the largest off- 
centering aside from Lu are of the La ions and depending 
on the supercell these may or may not be parallel to the 
Lu. 

The individual Lu off-centerings tend to avoid [111] 
and equivalent directions, and also tend to be non- 
coUinear with each other for 50% Lu concentration, e.g. 
in the "G" ordering (nominally the highest symmetry 
case), three of the Lu displace along different Gartesian 
directions, while the fourth has a smaller displacement 
near [111]. A preference for Gartesian directions was 
noted in (K,Li)Nb03)^ and understood from the fact 
that the square faces of the cage (the faces with the most 
room for the Li ion) are along these directions. Turning 
to the question of polar behavior, both of the cells at less 
than 50% Lu concentration show polar structures. At 
50% the cells "G", "P" and "A" have polar structures, 
most strongly so for "G", while "C" is nearly antifer- 
rodistortive, with tilts that avoid compressing the O - 
La bonds. Of the four supercells investigated the relaxed 
"G" structure had the highest energy. Taking this energy 
as the zero, the calulated energies were -0.43 eV, -0.36 eV, 
and -0.42 eV, for "A" , "C" and "P" , repectively on a per 
formula unit (10 atom) basis. The similarity of the en- 
ergies other than "G" show that the A-site ions will be 
disordered in material made by conventional methods, 
in agreement with what would normally be expected in 
perovskites with chemically similar, same charge, ions. 
While in the ordered "C" structure structure a long range 
tilt pattern of this type is allowed, this will not be the 
case in general for disordered A-site ions. 
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Assuming that the ^-sites are disordered in the al- 
loy, this dependence of polar behavior on chemical order- 
ing is more consistent with relaxor ferroelectric behav- 
ior than normal fcrroelectricityi^SiiSiiil However, it should 
be kept in mind that the relaxations reported here were 
performed within the LDA, which also predicts zero or 
very small band gaps, while in reality larger, but un- 
known gaps may be present. Larger gaps would lower 
the electronic dielectric constant favoring ferroelectric- 
ity and stronger coupling between the Lu off-centerings. 
While actual ferroelectricity may occur, what can be con- 
cluded here is that polar behavior will occur in disordered 
(La,Lu)MnNi06 for Lu concentrations at or below 50%. 
This may be ferroelectricity or relaxor ferroelectricity. 

This polar behavior arises because of frustration of the 
tilt instabilities due to the mixture of ^-site cation sizes 
and the fact that the coherence length for off-centering 
of j4-site ions is shorter than that for the tilt instabili- 
ties. This mechanism is quite general in principle, may be 
useful in producing polar behavior in other perovskites. 



Qualitatively, this is related to the rigidity of the BOq 
octahedra. This condition is often but not always met, as 
for example, while often tilt instabilities are strengthened 
by pressure, there are cases where this does not holdJ^ 
It may be difficult to synthesize perovskite 
(La,Lu)2MnNi06 due to phase separation^ or compet- 
ing phases, e.g. tungsten bronze as often occurs with 
mismatched A-sites. These issues can sometimes be 
overcome using thin film techniques, such as pulse laser 
deposition, or by high pressure synthesis, which favors 
the high density perovskite structure. Further, polar 
behavior is predicted over a wide composition range. 
This may help in finding specific compositions amenable 
to synthesis. In any case, the proposed mechanism is 
quite general, and should apply to other mixtures of 
A-sitc ions with different size, e.g. La with other small 
rare earths. 
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